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ABSTRACT 

To discuss the quantum to classical transition in quantum cosmology, we study 
the decoherence factor and the peak of the Wigner function, which respectively 
represent the degree of decoherence and the degree to which the classical motion 
of the Universe is defined, in a Robertson- Walker universe model coupled with a 
scalar field. It is known that the decoherence factor is divergent in some cases. 
This implies that perfect decoherence occurs, and classical correlation criterion 
fails. In this paper we discuss the divergence of the decoherence factor in some 
detail and obtain the constraints for decoherence factor to be convergent, making 
use of the arbitrariness defining the reduced density matrix. The result is discussed 
in connection with the arbitrariness of system/ environment splitting. 



1. Introduction 



Many researchers accept that quantum theory can be applied to all events and 
even to the entire Universe. However, macroscopic objects, including our universe, 
behave classically. If we seriously consider quantum theory as a universal theory, it 
is a crucial problem to derive classical behavior from quantum theory under some 
suitable conditions. 

To characterize classicality, two criteria are generally used as fundamental, 
classical correlation and quantum decoherence. 

Strong correlation is necessary to predict from the wave function of the Uni- 
verse. It is a reason that we cannot have probabilistic predictions from the wave 
function of the Universe that is a single system and can only predict the event 
whose (conditional) probability is almost unity, i.e., the strong peak of the wave 
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function or Wigner function. ' The existence of correlation between configura- 
tions and momenta along a classical trajectory implies that the system follows 
classical equations of motion. In order to analyze whether a given wave function 
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has classical correlation, we often examine the peak of the Wigner function. 

We may expect that the WKB state is classical in a sense and that the Wigner 
function associated with it shows classical correlation. However, the Wigner func- 
tion of the WKB state not only has no peak around the classical trajectory in 
phase space but also oscillates very rapidly.' 4 ' 51 The cause of these phenomena is 
the existence of quantum interference. To acquire classical correlation, we need 

. . [6] 

some coarse graining. 

Quantum decoherence induces the coarse graining necessary for the Wigner 
function to have a peak! 61 And originally, quantum decoherence is necessary for in- 

ryi 

dividual systems in a composite system to have definite characters, and necessary 

rgi 

to assign the probability into each history of the system. 

Classical correlation and quantum decoherence are not only complementary but 
also exclusive in part, because decoherence induces a coarse graining on the Wigner 
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function, and the coarse graining leads to a spreading of the Wigner function. 
Let us assume that we have a superposed state in the position representation, 
and after the decoherence mechanism acts, its density matrix associated with the 
state is diagonalized in that representation. Therefore we have a classical ensem- 
ble consisting of the states whose positions are decided. Each state is a position 
eigenstate so that the conjugate momentum is indefinite. So we have no classical 
correlation between positions and momenta in phase space. To obtain a definite 
classical correlation, decoherence should be moderate. 

A mechanism leading to decoherence is the "environment induced superselec- 
tion rule"! 101 Macroscopic objects are not isolated from their environment. We 
finally need reduced information because we cannot observe the entire Universe. 
This inevitable interaction with the environment induces non-unitary evolution on 
the relevant degrees of freedom called the system. 

This mechanism has been applied to quantum cosmology! 11 161 There, decoher- 
ence of the spacetime and classical correlation, i.e., derivation of the semiclassical 
Einstein equation, has been discussed. Decoherence and correlation of the scale 
factor of the Robertson- Walker universe by an inhomogeneous scalar field as envi- 
ronment are usually examined under the WKB ansatz, and the degree of decoher- 
ence is expressed by the overlap integral, called the decoherence factor, between 
the states of scalar fields which have evolved on different WKB histories. 

In some cosmological models, we face perfect decoherence of the scale fac- 
tor! 121151161 Because we deal with a system of an infinite numbers of degrees of free- 
dom, we essentially need regularization and renormalization procedures. However, 
this feature is not removed by the standard renormalization procedure! 161 Further 
this is an undesirable feature for classical correlation. So we usually introduce a 
cutoff to the degrees of freedom of scalar modes as the environment. Unfortunately, 
whether the scale factor behaves classically is strongly dependent on the choice of 
cutoff. [nI12I61 Meanwhile, we have a model that shows moderate decoherence. Several 
authors related the cause of perfect decoherence to the lack of adiabaticity of the 
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WKB histories. 1 1 J 

[17] 

We always have arbitrariness in defining the reduced density matrix. The 
reduced density matrix may be defined by integrating out the modes {4> n } as the 
environment, 



Alternatively, we may define it by integrating out {f n — dg( a ) ( l ) n} as the environ- 
ment, 



The arbitrariness of the choice of the environment variables brings change to the 
conjugate momenta so that the portion spanned by the Hilbert space of the system 
within the total Hilbert space changes. Except for the coincidence limit, a = a', 
Pred 7^ Pred- Then, is it a natural choice which reduced density matrix we use for 
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the decoherence argument of the scale factor? Laflamme and Louko assert that 
the choice g(a) = 1 is preferable to g(a) = a^ 1 because the decoherence induced 
by the former is more directly related to particle production. 

In this paper, through a model consisting of a scale factor and free scalar 
field with curvature coupling, we somewhat carefully discuss the question of which 
which choice for reduced density matrix leads to moderate decoherence. The paper 
is organized as follows. In Sec. II we define the decoherence factor induced by inte- 
grating out the rescaled scalar field according to the usual procedure. In Sec. Ill we 
relate the variance of Gaussian solution for the Schrodinger equation of rescaled 
scalar field to mode functions that are solutions of Klein-Gordon equation, and 
review the known results briefly. In Sec. 4, using the constraints on large frequency 
behavior of the mode function by Bogoliubov implementability, we obtain the con- 
straints on rescaling g(a) for the scalar field and the curvature coupling constant. 
Sec.V is devoted to a summary. 
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2. Decoherence Functional 



We consider the Robertson- Walker universe with a free scalar field, and exam- 
ine the decoherence of the scale factor by using the rescaled field as the environ- 
ment. The Lagrangian density of this system is 



£ = 



M 2 
12 



-9 R 



-9 



[(V$) 2 + m 2 $ 2 + £R<S> 



(2-1) 



We choose rescaled modes {f n } as the environment, 



„ a 9W 



(2.2) 



where {Y n } are normalized eigenf unctions of the Laplacian with respect to hij, and 
the index n is an abstract index classifying the eigenf unctions Y n . The eigenvalue 
of Y n is -n 2 = — (k 2 — K). 

Using the RW metric, 



ds 1 = -N{tydt l + a\t)hij{x)dx l dx 3 

= - iv2 * 2 +« 2 (i^ + r2 < iSi2 )- 



(2.3) 



the Lagrangian is written as 



G 



1 

2N 



—G^X^X" - NU , 



X» = (a,{fn}) , 



U{X) -~ M — + 2a^ fn 



~,2 



M 2 V(a) + -^-J 2 , 



lag 2 ' 



G n ,a — 



-M 2 a , 



E( Ma\ ri 
n \ g 



'2i 



M 2 a 2 



A, 



Qaa _ 



1 



a A- 



M 2 a ' 



(2.4) 



(2.5) 
(2.6) 
(2.7) 
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A„=G"" (in Jr^Sn, (2.8) 



M 2 a 2 g 2 (a) 



n 



where G aa is the supermetric in the case of no coupling to the scalar field, and the 
dash represents the derivative with respect to the scale factor a. 

The Hamiltonian of this system becomes 

H = hot^P^Pu + U , (2.10) 

/ Q aa F FX n \ 

G» v = 2 . (2.11) 

Quantizing this system, we obtain the Wheeler-DeWitt equation, 

[—□G + = , (2.12) 

where is the D'Alembertian with respect to supermetric G^, and the integral 
measure for the square of wave functions is da U n df n \f—G ~ rfa^/— G~^,H n [df n (a 1 / 2 / g)\. | 
Due to the existence of the Planck scale, we look for the solutions under the Born- 
Oppenheimer approximation, 

*(a,{f n }) = C(a)e iS W X ({f n };a), 
0(S(a)) ~ 0(M 2 ) . 

Substituting Eq.(2.13) into Eq.(2.12) and expanding in powers of M -2 , we obtain 
the Hamilton- Jacobi equation from leading term, 

]-G aa {d a S) 2 + M 2 V = , (2.14) 

To next order in M -2 , we obtain the probability conservation equation and the 
Schrodinger equation, 

d a (V^GG aa C 2 d a S) = , (2.15) 



6 



1 a 1 / 2 
t(G aa d a S)(d a + -d a \n ) X 

^ 9 



1 ~ 2 



2a n ' 2 
6e(l-60(G aa ^) 2 , 



(A n <9n + d n \ n ) + " 2 / n 



2a^ 



(2.16) 
(2.17) 



where we choose the prefactor C to be WKB prefactor in the case of the pure RW 
universe. By this choice of C, wave function of the scalar field is normalized as 



,1/2 



We define WKB time as 



/ n{<rK 



_L |. = (G aa <9 a S)<9 a , 
Ndt v ; 



(2.18) 



so G aa d a S = a/N. Eq.(2.16) is rewritten as 



1 dv> 

N~dt 



a 2 1 

2a n + 2ag 2 nJn 



at 



2 a ( ag\ ( agy 



(2.19) 
(2.20) 

(2.21) 



In quantum mechanics, superposed states can also exist, and the general form 
of wave function is 



* = ^2Ci(a) exp[iSi(a)}xi({f n } ; a) , 



(2.22) 



where the index / is a parameter that classifies the solutions of Eq.(2.14) . In 
our case, I = ± where +(— ) corresponds to expanding(contracting) universe, 
S- = —S + , and d/Ndt- = —d/Ndt + . Hereafter we use a "dot" to represent 
the derivative with respect to t+. 
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We define the density matrix of the total system and the reduced density matrix 
of the scale factor. First, we define the total density matrix as 

p(a, {f n } ; a', {f n } ) = *(a, {f n } {f' n } ) . (2.23) 

From the bi-scalar character of the density matrix, we define the reduced density 
matrix as 

f r /o 1 / 2 \ 1/2 /a' 1 / 2 \ i/2i 
Pred(a,a') = j u{df n {^) ( ) }p(a, {/»} 5 a', {/»} ) ■ (2.24) 

From Eq.(2.22) , we obtain 

p red (a, a') = Pi°i> (°> a ') h,i> (a, a') , (2.25) 

p$(a,a') = Q(a)^(a')exp[^(a) - S,,(a'))] , (2.26) 
IlfM = J ^dfneM-^l-0'i')fn\M{fn} I a)# ({/«} I a') , (2.27) 

^ ) = i^( ln ^)'' (2 ' 28) 

where is the density matrix in the case of a pure RW universe, and / is often 
called the decoherence factor. From normalization, Eq.(2.18) , at the coincidence 
limit, I = I' and a = a', we obtain \Iu(a, a)\ = 1. Further, by the Cauchy- 
Schwartz inequality we generally obtain the inequality |i///(a, a')\ < 1. Thus, the 
decoherence factor reflects the degree of decoherence of the scale factor. 
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3. Solution of the Schrodinger equation and Known Results 



As we discussed in the previous section, to examine decoherence under the 
WKB ansatz, we first solve the Hamilton- Jacobi equation, Eq.(2.14) , for the 
background geometry, and next solve the Schrodinger equation, Eq.(2.19) , in the 
background geometry that is a solution of the Hamilton-Jacobi equation. Finally, 
we estimate the decoherence factor, Eq.(2.27) . 

Because there is no mode coupling in Eq.(2.19) , we deal with each mode 
separately. Eq.(2.19) has a quadratic form, and we seek solutions of the Gaussian 
form 

1>rd = (^) 1/4 exp[^(t) - \B nl {t)fl\ , (3.1) 
where B n iR = Re(B n i) > 0. From the Schrodinger equation Eq.(2.19) , we obtain 

if = ~f a B M , (3.2) 
We introduce the variable u n \ as 



Ng 

and Eq.(3.3) is rewritten as 



Bnl = ^TT^l^Kl + ln #]' > ( 3 - 4 ) 



N\N Unl ) +u nUnl = U , 

vl = k 2 + a 2 [m 2 + (Z-±)R] . 



(3.5) 



Eq.(3.5) is just the Klein-Gordon equation. So u n i is a mode function, and the 
choice of the mode function corresponds to the initial condition of the wave func- 
tion, ^) n i. 
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In the Gaussian form, Eq.(3.1) , the decoherence factor, Eq.(2.27) is written 

as 



Iij>(a,a) = Uexp[i(A n i(a) - A nV {a'))\ 



^B n i R (a)B n i IR (a') 
(B nl (a) + B* nl> {a')f 



1/4 



(3.6) 



B nl (a) = B nl (a) + i 0,(a) = -i-^[\nu* nl + (6£ - 1) In a]' . (3.7) 



Let us try to apply the procedures described above to some models. First, we 
consider the case of a massless minimally coupled field on the deSitter background 
as a zeroth order WKB solution! 11 ' If we choose the deSitter invariant vacuum as 
the boundary condition of the mode functions, the functions B n \ are given by 



B nl (a) = (k - l)a fc2 + R 2 a 2 _ 1 . 



(3.8) 



where we chose the rescaling g(a) = a ^ So we obtain 



"I 



\Ill,(a,a!)\ = n 



i (a 2 - a' 2 ) 2 (a' 2 WH' 2 a 2 - 1 - a 2 lWH' 2 a'' 2 - 1 ) 2 i 
Aa 2 a' 2 Ak 2 a 2 a' 2 



(3.9) 

By calculating the infinite product, we obtain perfect decoherence between the 
expanding universe and contracting one, and perfect decoherence within a WKB 
branch, 



\I+-(a,a)\ = n 

k=2 

oo 

\Il,l{a,a')\ < n 

k=2 



1 + 
1 + 



2„2 



k 2 



-fe 2 /4 



a 2_ a /2x2 fc2 /4 



4a z a 



2 n l2 



= 0, 

= for a 7^ a' . 



(3.10) 



Next we consider the case of a massive conformally coupled field on an adi- 
abatically evolving universe! 10 " 161181 Let us set, at scale factor a = ao, the initial 
condition of the scalar field as in- vacuum and choose the rescaling g(a) = 1 and 
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N = a. Using Bogoliubov coefficients, which are estimated by a standard WKB 
technique, the in-vacuum is related to the adiabatic out-vacuum as 



u out ~ 

""n 



uj? = a n u° n ut + (3 n u° n ut * , 
exp(-i f' ' u n dr)') 



2 



Pn = ^ exp[-^(/c 2 + m 2 al)} , 



Otn ~ 1 ■ 
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7T 



\Zm 2 ao^'(ao) ' 



(3.11) 



Thus 



5 n+ (a) = i/„ 



out* , 



* /9* / 0'< 

a n ~ Pn u nl u n 
ru* -X- R*o,out l ,out* 



+ i 



(3.12) 



The norm of the decoherence factor between the expanding universe and contract- 



ing universe is given by 



[18] [16] 



OO 

\I+-(a,a)\ = n 

fc=2 



l-^{2|/5 n |cos(2 / Vn)-^\ ] 



2-i fc 2 /4 



(3.13) 



and due to 0{y n 2 ) = 0(z> n /V 2 ) = 0(k 2 ), the infinite product on the R.H.S. of 
Eq.(3.13) is convergent, so that moderate decoherence results. Similarly the deco- 



herence factor within a WKB branch is also calculated, and it does not vanish. 



[15] 



4. Constraints on rescaling and curvature coupling constant 

We would like to know how the scale factor decoheres. Thus we concentrate 
on the norm of the decoherence factor, 



\Il,l'(a,a) 



n 

n 



n 

n 



^B n i R (a)B nVR {a') 
\B nl (a) + B*(a>)\ 2 



1/4 



'M2-, 



1 - 



\B nl (a) + B* nll (a')\i 



1/4 



(4.1) 



By virtue of a well known theorem for infinite products, the divergence of Eq. (4.1) is| 
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equivalent to that of 



v |£Utt)-^'K)l 2 

^ \B nl (a) + B* nl ,(a')\* 



(4.2) 



In our case, in the large k limit we can replace the summantion for mode k into 
an integral over k: 

,2 



E 



dkk z 



If Eq.(4.2) is convergent, it must hold that for large k, 

\B nl (a)-B nl ,(a')\ 2 



\B nl (a) + B* nl ,(a 



AI2 



o(k~ 



(4.3) 



In order to obtain the constraints for the decoherence factor to be finite, we 
attach some physically reasonable conditions to the mode functions. We expect 
that at each time we can have a particle picture and that the Fock representation set 
at each time can be related unitary equivalently, i. e., Bogoliubov implementability. 

This Bogoliubov implementability gives some constraints on the behavior of 



the mode function in the large k limit: 



[20] 



Unl = Re 



—i 



v n + oik' 



In 



-Im 



. la u* ' 



Nu 



* 

nl 



oik- 1 / 2 ) . 



(4.4) 

(4.5) 
(4.6) 

(4.7) 



For example, the behavior of each Eqs.(4.4) and (4.6) in the large k limit for 
adiabatic vacuume is, respectively, 



linl = v n + o(k 2 ) , 
Inl = o(k~ 2 ) , 



(4.8) 



so that Eqs.(4.5) and (4.7) hold. For the Hamiltonian diagonalization vacuum, 
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the behavior of Eq.(4.6) in the large k limit is 



la 

lul = (1-60^, 



(4.9) 



, [19120] 



so that Eq.(4.7) does not hold, except for when £ = 1/6. 

Under the WKB ansatz, Eq.(2.22) , we take the constraints Eqs.(4.5) and 
(4.7) for granted because the back reaction correction to the background evolution 
should be a small finite quantity so that the particle production is finite. 

Due to the constraints (4.5) and (4.7) , 



B 



nl 9 

g 2 



(4.10) 



and thus, from Eq.(4.3) , we obtain the constraint 



\B nl {a) - B nV {a')\ 2 ~ oik' 1 ) . 



(4.11) 



From Eq.(4.11) , we derive the constraints on rescaling g(a) and the curvature 
coupling constant to obtain moderate decoherence. 

4.1. DECOHERENCE BETWEEN DIFFERENT WKB BRANCHES 



As different WKB branches in our model, we have the expanding universe and 
contracting universe. First, we examine the decoherence between these two states. 
Its degree is measured by 



|i+-(a, a) | = n 



1 - 



\B n+ (a) - B n -(a)\ 



1/4 



(4.12) 



\B n+ (a) + B*_(a)|2 
and its convergence is equivalent to Eq.(4.11) with a' = a: 

\B n+ {a)-B n -{a)\ 2 = \ [( /Un+ -/i n _) 2 + { 7n+ -7 n _-2^(lna 1 -^)'} 2 ] ~ o{k~ x ) . 

9 (4-13) 
Under Eqs.(4.5)and (4.7) , for the decoherence factor, Eq.(4.12) , to be finite, we 
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obtain the constraint on the curvature coupling constant 



e = \, (4.i4) 

i.e., conformal coupling. 

4.2. DECOHERENCE WITHIN A WKB BRANCH 

Next we examine decoherence of the scale factor within a WKB branch. We 
omit the index / in the decoherence factor, and so on. The degree of decoherence 
is measured by 

|/(a,a)l = n[l- |Bn(o)+B , (o , )|2 j . (4.15) 
Its convergence is equivalent to Eq.(4.11) with I' — I, and 



|B»-^(«')| 2 =(^-^) 2 
\g 2 (a) g\a')) 



+ 



7n(a) ~ fQna 1 ^)- JnW) ~ irilna' 1 -^ )-^ (416) 



Under Eqs.(4.5)and (4.7) , for the decoherence factor, Eq.(4.15) , to be finite, we 
obtain the constraint for rescaling, g(a), and the curvature coupling constant: 

g(a) = constant , 

1 ,4,7) 
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5. Summary and Discussion 



In the previous section we obtained the constraints on the rescaling g(a) and 
the curvature coupling constant to make decoherence moderate irrespective of the 
evolution of the scale factor. Therefore, if the constraints (4.17) do not hold, we 
cannot obtain moderate decoherence, even with adiabaticity of the evolution of 
the scale factor, except for in a static region. Even though we choose the initial 
condition of the state of the scalar field to be a physically reasonable state leading 
to finite particle production, we have perfect decoherence unless the constraints 
(4.17) hold. So the divergence of the decoherence factor is not attributed to the 
lack of adiabaticity of the zeroth order WKB solution, but to system/environment 
coupling and system/environment splitting. The former is of dynamical origin, but 
the latter is simply of kinematical origin. 

Our result, Eq.(4.17) , is understood as follows. Because we deal with an 
environment with an infinite number of degrees of freedom, for the decoherence 
factor not to vanish, both xi( a ) an d Xl'( a ') must belong to the same Hilbert space. 
In other words, the states of almost all modes in xi( a ) should be nearly equal to 
the corresponding modes in xi'( a ')- Let us consider the case that \Iiji(a,a')\ = 1, 
so that the scalar field essentially has no coupling to the scale factor. This is 
the case when we choose g(a) = 1 for the massless conformal coupled scalar. 
Unless we choose g(a) = 1, the effective frequency of each mode is Q e ff = k/g (a) 
and the scale factor strongly couples to modes irrespective of k, d a fl e ff/fl e ff ~ 
0(k°). Even if we choose g(a) = 1, in the case of non-conformal coupling, we have 
momentum-momentum coupling, and this gives not only a-dependence but also 
direct /-dependence through d a S[ to the covariance, B n i(a), of the wave functions 
of all modes. Therefore, the choice of g(a) ^ 1 and non-conformal coupling induce 
certain differences to all modes, so that the overlap integral vanishes. Meanwhile, 
the mass term is ineffective for large k modes and it has no constraint. 

We calculate the quantities among a wider class in discussing decoherence 
than in ordinally discussing field theory on a curved space. On the decoherence 
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argument, we calculate the overlap integral between the states with the same initial 
condition at different times and in different universes. On the other hand, for 
example, particle creation is calculated by making a comparison between different 
states at the same time and in the same universe. In usual arguments for field 
theory on a curved space, rescaling g(a) and £ are ineffective. For example, in 
discussing that two particle models belong to the same Hilbert space, we, in the 
Schrodinger picture, calculate the decoherence functional for I' = /(same universe) 
and a' = a(same time), but for different particle models({w n } and {v n }): 



\I(a, a) 



n 

n 



\BJa) 



B' n (a)\ 2 



\B n (a) + B'*(a) 



1/4 



(5.1) 



where B n and B' are written as 



B n (a) = 
B' n {a) = 



la 



Ng 2 (a) 



la 



Ng 2 (a) 



ln<(a) +lna 1_6 ^ 



(5.2) 
(5.3) 



Therefore we obtain the well known result 



\I(a, a) 



n 

n 



n 

n 



\U*V* - wVlV/4 



\u*v — u*v\ 2 



\Oir 



1/4 



n|a n r 1/2 , 



(5.4) 



where v n = a n u n + fi n u* n and the index / is omitted. Thus by requiring Eq.(5.4) to 
be convergent, we obtain the constraints only on mode functions, and these are 
nothing more than Eqs.(4.5) and (4.7) . 

In an anisotropic universe, even if we choose g(a) = 1, modes tightly couple 
to geometry irrespective of k, dVt n /VL n = [d{e 2 ^)ijk l k^ + ..\/[{e 2 ^)ijk l k^ + ...] ~ 
O(k ). We believe that in anisotropic case, the decoherence factor must vanish. 
Further, it is known that there is no particle model which leads to finite particle 
production.' 191 This may make the WKB ansatz, Eq.(2.13) , invalid. 
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As described in the Introduction, we have the arbitrariness of defining the 
reduced density matrix, i.e., one of system/ environment splitting. At least in 
isotropic case, the arbitrariness is important to obtain the moderate decoherence 
that is necessary for the classical correlation criterion to hold. By requiring mod- 
erate decoherence, we can obtain a system/environment splitting, g(a) = constant, 
which is characterized by no coupling between the geometry and scalar field in the 
large k limit. 

If we really need moderate decoherence, it may fix a part of the arbitrariness 
of the system/environment splitting in a sense. 
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